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PsJ ■ Abstract 

j^ I Understanding the relationship which integrable (solvable) models, all of 

which possess very special symmetry properties, have with the generic non- 
integrable models that are used to describe real experiments, which do not have 
C^^ ■ the symmetry properties, is one of the most fundamental open questions in 

both statistical mechanics and quantum field theory. The importance of the 
(~| ■ two-dimensional Ising model in a magnetic field is that it is the simplest system 

CIhI where this relationship may be concretely studied. We here review the advances 

made in this study, and concentrate on the magnetic susceptibility which has 
revealed an unexpected natural boundary phenomenon. When this is combined 
with the Fermionic representations of conformal characters, it is suggested that 
the scaling theory, which smoothly connects the lattice with the correlation 
length scale, may be incomplete for H ^Q. 
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It may be rightly said that the two dimensional Ising model for iJ = is one 
of the most important systems studied in theoretical physics. It is the first sta- 
f~j \ tistical mechanical system which can be exactly solved which exhibits a phase 

PsJ ■ transition. From the exact results for the free energy [1], spontaneous magneti- 

zation [2], [3] and correlation functions [l]-[8] a point of view has been developed, 
which embraces the concepts of scaling, universality and conformal field theory, 
that extends the exact results of the Ising model to more general situations. 
These concepts are widely used to analyze both experiments and models of crit- 
j^ , ical phenomena. Furthermore the correlation functions provide very concrete 

realizations of the concepts of mass and wave function renormalization used to 
define Euclidean quantum field theories. 

However, starting with the work of Nickel [9], [10] on the magnetic suscep- 
tibility new properties of the Ising model have been uncovered |ll)-|27j which 
go beyond what has been seen in the computations of the free energy, spon- 
taneous magnetization and correlation functions. These new features need to 
be explored to see if there is relevant physics which is not incorporated in our 
current view of critical phenomena. In this article we will review these new phe- 
nomena and the relation they have with scaling theory and Euclidean quantum 
field theory. 

In section 2 we define what will be meant by an Ising model. In section 3 



we review the known exact results ior H — 0. In addition to the well known 
results for the free energy [1] and the magnetization [2], [3] we will put partic- 
ular emphasis on the magnetic susceptibility which has an expansion analogous 
[5],[5],[Tn] to a Feynman diagram expansion. These Ising model integrals share 
with Feynman diagram integrals the property that the integrands are algebraic 
functions. It was shown, long ago, by Kashiwara and Kawai [25], that these 
Feynman integrals are holonomic (i.e. they satisfy overdetermined systems of 
partial differential equations). However, the infinite sum of diagrams will not 
have this propertjQ- These problems of power series expansions, where the in- 
finite sum has different analytic properties than the coefficients of any finite 
power of the expansion parameter, are seen in the magnetic susceptibility of the 
Ising model [29]. When expanded about T = Tc, this expansion of the suscepti- 
bility is not convergent but is only asymptotic 11 [26 . This feature also occurs 
in Quantum Electrodynamics |30) and Quantum Chromodynamics [31] . 

In section 4 we extend our considerations to iJ 7^ 0, where much less is 
known and there is much to be learned. We will present perturbative studies, 
for H near zero, of the two-point function [32], and the two-particle scattering 
amplitude [33] . We will also present integrable perturbations of conformal field 
theory [31]-[SS] about T = Tc,H ^0, about the Lee- Yang edge [37], [55], and 
the relation to Rogers- Ramanuj an identities [3S]-[13]. We will finally discuss 
scaling theory [TT] , [33] , [33] , [3S] , and the study of the general non integrable 
perturbation of [3S]- 

We conclude, in section 5, with an answer to the question of "Why is the 
Ising model is important?" 

2 What is the Ising model? 

We begin by defining what we mean by the two-dimensional Ising model in a 
magnetic field. 

The simplest, and most well known, case is for nearest neighbor interactions 
on a square lattice defined by the classical interaction energy 

£i{H) = - ^ ^ {E"" ■ aj^k cTj+iM + E^ ■ aj^k (Tj,k+i + H ■ aj,k}, (1) 

where Uj.k = ±1 specify the "spin" at row j and column A; of a square lattice. 
This definition can be extended to nearest neighbor interactions on other lattices 
by a straightforward expansion of the notation. We will impose either periodic, 
or cylindrical, boundary conditions. 

For this interaction energy we are interested in computing the partition 
function 

Zi„,i.(T,i/) = ^ cxp(-£:,(i/)/A:BT), (2) 

<T = ± 



^It has first remarked by Enting and Guttmann 1291 . on the anisotropic lattice, that the 
susceptibihty of the square Ising model is not D-finite 1291 . The isotropic susceptibility is also 
non-holonomic (but this is not a consequence of the previous anisotropic result). 



where /cb is the Bohzmann constant and the sum is over all values of the spins 

From the partition function ^ one gets the free energy F(T, H) in the 
thermodynamic limit 



F{T,H)/kBT = hm ^lnZi„,i.(T,iJ), (3) 



the magnetization 

the magnetic susceptibility 

x(r,ij) = ^^^^^ = i^ E E {KoaM.A^)->('(r,i?)}, (5) 

and the internal energy 



dH keT 



u = kBT^—F{T,H)/kBT, (6) 

where we have used the definition of the thermal average of an operator O 

(O) = lim Z-}^,{T,H)- y20cxp{-£i{H)/kBT). (7) 

We note, in particular, that the magnetic susceptibility, at _ff = 0, is written, 
in terms of the two-point function, as 

OO CXD 

where 

M = lim M(T,H), (9) 

is the spontaneous magnetization. For the nearest neighbor interaction ([T]), the 
internal energy ((6]) reads: 

u= -E-" ■ KoCTi.o) -^''- KoCTo,i) -^- Ko). (10) 

The interaction energy of nearest neighbor Ising model ([T]) may be general- 
ized to interactions farther than nearest neighbors with interaction energy: 

L" L'^ L" L** 

£g = - Y. T. E E i?(b-/Ufc-A:'i)-a,,fca,',fc, 
-^- J2 T. ^^.'=- (11) 

j^-L'^ k=-L>' 



This reduces to © when £^(1,0) = E"" /2, £'(0,1) = E^ /2 and all other 
E{j^ k) = 0. The free energy, magnetization and susceptibihty are all computed 
by replacing £j in ^ hy Eg- The internal energy pU)) generalizes to 

" = - E E ^(ijUfci) • Ko^,,fc) -ff-Ko). (12) 

3 What do we know for H = ? 

We have obtained more exact results for the nearest neighbor Ising model at 
H = than for any other system in statistical mechanics. 

3.1 Free energy 

The free energy of the nearest neighbor Ising model with H = was computed 
by Onsager [T] in 1944 

-F/keT = ln2 (13) 



2iT cItt 



-^ / d^i / d02ln 



, 2E^ , IE"" 2E^ ^ , 2E'" 

cosh — — ; cosh - — — - smh - — — cos Oi ~ smh — — ; cos O2 
ks-L Kb! kb-l ks-L 



This free energy has a singularity at T = Tc determined by 

sinh2£;''/fcBTc- sinh2£;''/fcBTc = 1. (14) 

This temperature Tc is referred to as the critical temperature. The free energy 
may be expanded, about T = T^ as 

-F/kBT ^ (T-r,)2. ln|r-r,|- i^i(T-T,) +F2{T-T,), (15) 

where £1, and £2, are analytic at T = Tc (i.e. they both have convergent power 
series expansions). 

3.2 Spontaneous magnetization 

The spontaneous magnetization was first announced by Onsager [2, in 1949, and 
a proof was given by Yang [3 in 1952 

M = (l-fc2)i/8^ (16) 

where 

yfc = (sinh2£;7fcBr- sinh2£;^/fci3r)"\ (17) 

for T < Tc, and is zero for T > Tc- The history of the Onsager result is given 
in a most interesting paper of Baxter [4 7) . 



3.3 Correlation functions 



All correlation functions of the nearest neighbor Ising model ai H = may be 
expressed as determinants. These are particularly simple for the row correlation 
(co.oCo.iv)? and the diagonal correlations (o'coO'a'.a'), which are given by 



Dn = 
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with 
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de 



a-N+i 
a-N+2 



(l-aie*^)(l-a2e-*«) 



-,1/2 



(18) 



(19) 



where for {(To.o<^n,n} 

ai =0, a2 = {smh2E'' /keT smh2E'' /keTy'^ , (20) 

and for (cro,ocro,jv) 

ai = e-^^^/''^^ ■ tanh E'^/keT, aa = e-^-EVfe^T . coth E'^/keT. (21) 

From this determinantal representation (IT51) of (cro,oO'o,Af)i and {cro,o<^N,N), 
it is easy to obtain the behavior of the correlation function as T — >■ Tj,. The 
integrals a„ all have logarithmic singularities at T = Tc- and we find [TT], [JS] 
for both the row and the diagonal correlations that Dn has an expansion of the 
form 



N oo 



(22) 



Dn = J2Y1 rf('''")(ln|r-r,|)P. \T-n\P+", 

where for each p the sum over n converges (and thus defines an analytic func- 
tion). The correlations {(Jo.o(Jm,n) have a similar expansion in terms of powers 
of In |T — Tc I multiplied by functions which are analytic at T — Tc. 

When T — Tc the determinant for the diagonal correlation reduces to a 
Cauchy determinant. It is explicitly evaluated as |5], yielding 



■2\^ 
[o'o,q<^n.n) = I — 



From this, we find as A^ - 

{o'o,oO'n,n) 
with 



CXI that 
Ac 



Ai/4 



N-l , ^ l-N 

■n('-4? 



o[N- 



(23) 



(24) 



(25) 



However, for T ^ Tc, the representation ([18)) is not an efficient way to study 
the limit M,N ^ oo. 



64A^ 



Ac = 2^/12 . exp[3C'(-l)] -- 0.6450 • ■ 
where C'i^) is the derivative of the zeta function. 



3.4 Form factor expansions 

To study the correlation functions {oQfiaM,N) when M,N — >■ cxi, the determi- 
nants are re-expressed in an "exponential form" 6_ for T <Tc &s 

oo 

KoaA/.Ar) = A^'- exp^F!")(r;M,7V), (26) 

n=l 

and, for T > T^ 

OO C« 

(ao,oaM,JV> = Ml-Y, G^^^-^\T-M,N) ■ exp^ ^'"'(T; A/,7V), (27) 

j=l n=l 

where 

M+ = [1 -(sinh2£;7fcsr- sinh2£;'Vfcs^)Y^^ (28) 

is referred to as the disorder parameter and is the value of the spontaneous 
magnetization on the dual lattice where E"" and E^ are replaced by E*^ and 
E*"" defined as 

sinh2£;*7fcBr = l/sinh2^7fcBr, smh2E*'' /ksT = 1/ smh2E^ /keT. 

(29) 
The functions F^'\T]M,N) and G^"-\T;M,N) are n-dimensional integrals. 
These exponentials can be expanded, and terms combined [6], [9], [10] into 
what is referred to as the form factor expansion. For T <Tc this expansion is 

oo 

('To.ocTA/,^) = (1 - tf" • {1 + E fM%{T)}, (30) 

n=l 

with fMiq{T) a 2?i-dimensional integral and: 

t^ {smh2E'"/kBT- smh2E''/kBTr^. (31) 

For T > Tc the form factor expansion reads 

oo 

{ao.o<7M,N) = {l-ty/'-J2&'^(T^^ (32) 

n=0 

where fj^/^j^ (T) is a 2n + 1 dimensional integral, and 

t= {smh2E" / keT ■ smh2E''/kBT f- (33) 

For the diagonal correlations {<Jo,o<^n,n} a simpler alternative form of fj^j^{t) 
is given in [T7] , and proven in [TH] . 

The behavior of (crOjOcrM.Ar) is easily obtained, from this form factor expan- 
sion when T 7^ Tc, because JmnC^) ^^^ ^^ exponential decay 

fl:^VT) - e--(^^*^/^)-«, (34) 



where i?^ = M^ + N'^, where k{T;M/N) depends on the ratio M/N, and 
vanishes, when T —^ Tc, as: 

k{T-M/N) ^ \T-T,\. (35) 

This exponential decay defines a second length scale which, when T — > Tc, is 
infinitely great compared to the lattice length scale which defines the interaction 
energy ([IJ. It is worth noting that the form (IMl) implies that the angular 
dependence of the decay of the correlation functions is the same for T > Tc and 
T<T,. 

However, as T — ^ Tc, each term of the form factor expansion vanishes because 
(1 — t)^/* vanishes, and each f^j^iT) is finite at T = Tc- The singularities of the 

/m Ar(-^) ^^^ ^^^ at T = Tc- These functions satisfy Fuchsian linear differential 
equations [T7], and have logarithmic singularities at T = Tc where the highest 
power of In \T — Tc\ is n. It therefore requires a detailed infinite summation of 
powers of logarithms to reproduce the behavior (|22|) for fixed M, N a.s T -^ Tc- 
In fact, this type of demonstration has never been carried out. 

3.5 Differential equations for (cro,oC"Af,Ar) 

The correlation functions (cro,o<''M,Af) ai'e holonomic (D- finite): they satisfy a 
set of partial linear differential equations in the variables sinh E^'^/ksT. This 
is exactly the holonomic property of Feynman integrals, shown by Kashiwara 
and Kawai [28]. For the diagonal case the more specialized result holds that 
(cro,oO'Ar,Ar) satisfies a linear Fuchsian equation. The order of this linear differ- 
ential equation is iV + 1, which is the minimum order needed to accommodate 
the singular terms at T = T^, In^ |T - T^ with 0<p<Noil^, which were 
directly obtained from the determinantal representation (J18p . A few examples 
are given in jl6) . 

However the diagonal correlation (o'coO'Af.A') has the much more remarkable 
property, discovered by Jimbo and Miwa [S] in 1980, that it satisfies a second 
order non-linear differential equation 



/ , , da 



For T > Tc the diagonal correlation is related to a by 

a{t) ^ t- (t-l)- -\og{ao,oaN,N) - ^, (37) 



with the boundary condition at i = 



(l/2)Ar 



(ao.o-Ar.A.) = t^/'- '-^ +0{t'+^/% (38) 



and, for T < Tc, it is related to a by 

a{t) = t- (t-l)- -\og{aQ,o<JN,N) -^, (39) 

with the boundary condition 

where {a)N = a(a+l) • • • (a+N—l) is Pochhammer's symbol {N > 1, (a)o = 1). 
These boundary conditions are obtained from the leading terms of the form 
factor expansions as i — >■ 0. Equation ((36)) is an alternative form of the Painleve 
VI equation [49] , called the cr- form of Painleve VI. 

No nonlinear isomonodromic (Garnier [50(l51j) systems of partial differential 
equations have been derived for the general correlation function {ao^oaM,N), 
even though they almost certainly exist. Such isomonodromic systems would 
yield, by differential algebra elimination, in the isotropic case, higher nonlinear 
order differential equations with the Painleve property. Some exact results for 
{o'o,o<^N,N-i) are actually given by Witte in |52) . 

However ((To.o'''J\/,Af) does satisfy quadratic difference equations [S3], [S3] 

sinh2£;''/fcBr • {C{M, Nf - C{M, N - 1) C{M.N + 1)} 

+ siTLh2E*''/kBT ■ {C*{M,Nf -C*{M -l,N)C*{M + l,N)} = 0, (41) 

sinh2£;7fcBr- {C{M,Nf -C(Af-l)C(M+l,iV)} 

+ sinh2^*'VfcijT- {C*(Af,A^)2 -C*(Af,A^- l)C*(M,iV + l)} = 0, (42) 

where 

C{M,N) - {aoflCTM^N), (43) 

and where C*{M,N) are the correlations on the dual lattice (|29]) . At T — Tc, 
where siTih2E*ykBT = sinh2EykBT ioT i = v,h, and C'*(M,iV) = C{M,N), 
these equations reduce to the discrete imaginary time Hirota equation |55] . 

3.6 Susceptibility 

By using the form factor expansions ([50]) and l\'32\i in ([S]) , we obtain the expansion 
for the susceptibility as the infinite sum of n "particle" contributions 

oo 

ksT ■ x+(T) = (1 - t)i/4 • t-V4 . J2 /''+'\T), for T > T„ (44) 

OO 

ksT ■ x-{T) = (1 - tY'' ■ Y^x^^^KT), for T < T,. (45) 

where explicit forms of the x (T), as j — 1 dimensional integrals for the general 
anisotropic lattice, are given in [11,. These series are convergent for T ^Tc- 



We may also consider what we call the "diagonal susceptibility", defined as: 

oo 

kbT-Xd{t) = ^ {KofTw,JV> -A^'}. (46) 



N=-cxi 



By use of the form of the diagonal form factor expansion of 53"[IH1 ^^ obtain 
as analogues of (^5]) and (03]) 

oo 

ksT ■ xMT) = (1 - 0'/' • t''^' ■ E xf +'^ W, for T > T,, (47) 

3=0 

OO 

ksT ■ Xd-iT) - (1 - tf'^ ■ E xf\t). for T < T„ (48) 

where the Xd (^) ^-""^ given in |51]. These series are also convergent for T =^Tc. 

The series expansions (HU and (|^5t are, what we referred to in the intro- 
duction as, the analogue of a "Feynman diagram expansion" . 

The integrals for x^"-'(T), and Xd (^)' have been extensively studied [H]- 
PS] . It is quite instructive to contrast them with the form factor integrals 
/m Ar(^) from whence they came. 

When T ^Tc each x^"^^) diverges as \T - T^-^ with a coefficient which 
depends on n and rapidly decreases as n increases. Thus, in terms of the variable 

T= -■ {amh~^2E/kBT ~amh2E/kBT), (49) 

we have, for the isotropic lattice [6], [56], [11] as r -> 

kBT-x±{r) ^ V2- C±-\t\-'/\ (50) 

The constants C_ and C^ are different, and are given as infinite series 

oo oo 

C_ = E ^*'"^' ^+ = E C''"+'\ (51) 

n— 1 n— 

where the C^") are C^") = 2"" 7r"-i !)„, with [57] 



2 



nlJo Jo wi Un (E"=i(%- + 1/"^))^' 

These integrals have been studied for n = 1, • • • , 6. The first terms in ([?T]) have 
been analytically evaluated in [6]: 

C(^' = 1, C(^) = -^. (53) 



The next leading term was analytically evaluated by Tracy [5S] as 

C(3^ = ^- (y +2 -3\/3Cl2(V3)) = 8.1446 •••X 10-^ (54) 
where 

aM ^ g!!^, (55) 

n— 1 

is Clausen's function. The next term read: 

C(4) = -^- (— -- --C(3)^ = 2.5448 ••• X IQ-^. (56) 

Accurate numerical evaluations have been made [57] for C*^^^ and C^^^ : 

C(5) = -^ X 0.0024846057 • • • = 7.9709118 • • • x 10"^ (57) 

C^^) = -^r-^ X 0.0004891422 • • • = 2.497501 • • • x 10~®. (58) 

A curious feature of these results is that the ratio C+/C- is found to be closely 
approximated by 127r, and the succeeding terms decrease by approximately three 
orders of magnitude. The study of the constants C_ and C+ has been contin- 
ued by high precision numerical computations [11', and the most recent evalu- 
ation [55], in 2011, is to 104 places. 

For the diagonal susceptibility each Xd (t) diverges as (1 — t)^^. One 
finds [27], [59] 

/I 2^l/4 "^^ 

kBT-x.+ {T) -^ ^'-y -E^lr^ (59) 

n—0 

where C^_^ — 1, and 

^'^+ = 3 ~ 9r2(5/6)r2(2/3) ~ r2(l/6)r2(l/3) " 0-016329 • • • (60) 
and 



ksT-xd- -^ {i-ty'/'-J^ct^^ (61) 



n=l 



where C^J = j, and 



8 I 37r2 



^^ = «■ l-^[64-16/]), (62) 



with: 

/ = -2.2128121 ••• (63) 



10 



In contrast to the form factors fj^j^{T), whose only singular point is T = Tc 
where Tc is real, the x (2^)'s have many further singularities. The first set of 
these singularities was found, by Nickel [9], [10], to be, for the isotropic case 
£■" = E'' = E, located at 

cosh^ 2E/kBT - smh2E/kBT ■ {cos{2Trj/n) + cos{2TTl/n) = 0, (64) 

with {[x] is the integer part of x) < j, I < [n/2], j — I — excluded (for 
nevenj + l~ n/2 is also excluded). Equivalentlv ([64|) reads sinh2£'/fcBTj^; = 
Sji = 1/2- (cos(27rj7n)+cos(27rZ/n)) ±j/2- [(4-(cos(27rj7n) +cos(27rZ/n))2]i/2. 
These Nickel's singularities are clearly on the unit circle \s\ — 1, or |fc| = 1. Do 
note that this is no longer the case for the anisotropic model. In the anisotropic 
case Nickel's singularities for the anisotropic x^^^'s become (see (3.22) of |11)): 

cosh^ 2E''/kBT ■ cosh^ 2E''/kBT (65) 

-(sinh2£'"/fci3r- cos(27rj/ri) +sinh2£;''/fcBT- cos(27r7n)) = 0, 

with j, I = 1, 2, • • • n. These (complex) algebraic curves (|65p . in the two 
complex variables sinh2£'''/fcBT, smh2E^ /ksT, are actually singular loco for 
the D-finite system of PDE's satisfied by the anisotropic x^")'s. These alge- 
braic curves accumulates with increasing values of n, in the same way Nickel's 
singularities ()64|) accumulate on the unit circle \s\ = 1, in a certain (real) sub- 
manifold S of the two complex variables smh2E'^ /kBT, sinh2E^ /kBT (four 
real variables). However, this "singularity manifold" S is not a codimension- 
one (real) submanifold (like the unit circle \s\ = 1 in the s-complex plane), 
but actually a codimension zero submanifold, as can also be seen on various 
analyses of complex temperature zeros (see for instance [60]). Note that this 
"singularity manifold" becomes very "slim" near the (critical) algebraic curve 
k = sinh2£'''/fcBT- smh2E^ /kBT = 1 (see for instance figures 1, 2 and 3 near 
the real axis in 60J). The relation between the analyticity in the two complex 
variables sinh2£'''/fcBT and smh2E^ /ksT, and the single complex variable T 
with E"" and E^ real and fixed, will be discussed elsewhere. 

Back to the isotropic model, introducing the variable 

"^ " 2- {smh2E/kBT +{smh2E/kBT)-^)' ^^^^ 

the singularities ([M)) for n — 3, 4, 5, 6 are given in table 1, where we note, when 
w is real, that sinh2£'/fcBT is real for —1/4 < w < 1/4, and is complex with 
I smh2E /kBT\ ~ 1 for \w\ > 1/4. Following [10] we define e, the deviation from 
the singular temperatures T-" determined by ([M|) . as s~^ — (1 — e) • sj^, for 
T < Tc, and s = (1 — e) • Sj,i, for T > Tc- Then, for T < Tc, the singularity 



in;(;(2")(T), at Tj,i, reads 



^(2«) . ^2„^-3/2^ (g7) 



■^This result on the anisotropic Ising model has been obtained from a Landau singularity 
analysis, generalizing I19II20| (S. Boukraa, S. Hassani and J-M. Maillard, unpublished results). 
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and, for T > T^, the singularity in x^^"+^H^) reads 



,(2n+l) 2n(n+l)-l 
^3,1 ' ^ 



Ine 



(68) 



The amphtude A", ; is given by [TO] 



with 



K 



(AT) 

J,' 



(iVzsin^,-0(^'-3)/' 



(sin2 0(i) cos(?!.« +sin2 0« cos<^(j'))(^'-i)/2 

rAf-l 



n^iK^^Va™) 



^(JV-3)/2)y]vr((iV2-l)/2)' 



^^^^ = ^-KrajN, 2cos6'j-, = cos (^^^'^ + cos (/)('\ 



(69) 

(70) 



Table 1: The Nickel singularities (|64| for the isotropic case E^ — E^ of x^"^ for 
n = 3,4,5,6 



n 



w 



-1/2, 1 
±1/2 



-l±v^ 



-1, 

±1, ±i/3 



3±75 
2 



The Fuchsian linear differential equations [12]- [IS] which the (isotropic) 
■^(") (7") satisfy, have singularities which accumulate on the unit circle, but, also, 
inside, and outside, this unit circle | sinh2_E/fcBT| = 1. However, the series ex- 
pansion of these integrals x'"^^) are actually analytic inside the unit circle 
(I sinh2i?/fcBT| < 1). The singularities of a linear ODE and the singularities 
of a particular series-solution of an ODE do not coincide. 

The singularities of the Fuchsian ODE may be obtained, from the integrand 
of the integrals, by the same "Landau" analysis [19] used to obtain singularities 
of Feynman diagrams. 

For the isotropic x^'^\ in addition to the unit circle singularities ([M]) . there 
are additional singularities at 



-3± i\/7 
8 ' 



where 



-1 ± i^/7 



(71) 



and they actually correspond to complex multiplication of elliptic curves, Heeg- 
ner numbers, and complex fixed points of the Landen transformation [20]. 

The singularities of x'"' (w) are to be contrasted with the diagonal suscep- 
tibility Xd (t) which only have [19] singularities on the unit circle \t\ = 1. 



12 



For T < Tc the singularities of xk " (^) ^''^'^ ^^ ^o = e^'^*-'/", and are of the 
form: 

A^^f ■ e2"'-i . Ine. (72) 

For T > Tc the singularities in Xd"^^\t) are at to = g^'^^/^^+Vz)^ and are of 
the form 

where the amplitudes A^^-' have yet to be determined. 

The linear differential operator for x^'^^ rightdivides the one for x*-^-* (in a 
direct sum structure [21p. Consequently, all the singularities of the linear dif- 
ferential operator for x*-^-* are also singularities of the operator for x'^''- It was, 
however, seen f21,, by means of a Fast Fourier Transform analysis of the series 
expansion of the integral for x S as a power series in t, that some x singulari- 
ties are not present in x^^^ • This is to be contrasted with a very recent result [HI] 
for the diagonal susceptibility x^ which has found that the singularities of Xd 
are present in xjj • 

3.7 Natural boundary for the isotropic model 

It is striking that the number of singularities (|64p in x (T) increases with n, 
and becomes dense in the limit n — ^ oo. This feature led Nickel [3] to the 
conclusion that, unless cancellations occur, there will be a unit circle natural 
boundary in the susceptibility x{T) in the complex s plane. The fact that ref.[2T] 
demonstrates that some singularities of the x*''^^ series are not singularities of 
the x''^'' series, supports the no cancellation assumption. The existence of a 
natural boundary, in the complex temperature plane, is not contemplated in 
the scaling theory of critical phenomena. The literature on natural boundary is 
quite narrow [62], as well as the methods and tools to study the neighborhood 
of a natural boundary. Curious situations may occur, like, for instance, the 
"radial porosity" one encounters with Borel's examples of monogenic functions 
(see section 10.5 page 21 of j63,). 

There is much that is, as yet, not understood about the properties of this 
Ising natural boundary. In the vicinity of any point on |s| = 1, such that s ^ ±1, 
the local spacing of the singularities, in x'""'\ is of the order n~^. However, for 
s = ±1 the local spacing is only n~^, and the dependence of the amplitude (|69|) . 
for large N, is different in these two cases. 

An initial analysis was made in [llj . for s — >■ 1, based on keeping only the 
singularities closest to r = 0. This analysis yielded an essential singularity 
of the form e~ '"^ where C is a constant. However, there is more to a nat- 
ural boundary than just an essential singularity, and further analysis will be 
required to fully assess the properties following from the singularities in x*^"^- 
In particular we note that, in general, the limiting value at a point on a natural 
boundary depends on the direction of the approach, and, as suggested in [11], 
an asymptotic expansion about |s| = 1, s ^ ±1, may not exist. The existence 
of a unit circle natural boundary in the isotropic square Ising model, seriously 
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questions most of the scaling arguments taken for granted on this model. For 
example it is desirable to explore the relation between the equality 7+ — 7_ 
for the susceptibility exponent 7, above and below Tc, and the unit circle nat- 
ural boundary which disconnect^ the inside and the outside of the unit circle, 
\k\ < 1 and |fc| > 1. 

3.8 Series expansions 

There is a second way of studying the susceptibility which is distinct from (and 
in a way complementary to) the use of the form factor expansion. By use of 
either (in principle) the determinant representation or (in practice) by the more 
sophisticated approach of using the nonlinear difference equations dH]), (|1^ 
with the explicit result for ((To,oO'jv,jv) of (P5|) . and for (cro,oCAr,Ar-i) of [S^, very 
long expansions about T — 0, and T = 00, can be obtained. For the isotropic 
lattice E"" — E^ = E this has been done in [11] and [26], and series expansions 
with may hundreds of terms have been obtained, not only on the square, but 
also on the triangular and honeycomb lattices. 

These long expansions have been analyzed in [TTj and [55] ■ This analysis 
leads to the following expansion of the susceptibility as T —>■ Tc which general- 
izes the leading diverging term given in ([50)1 



where 



and 



kBT-x± = C± ■ |r|-7/4 . F±(r) + B{t), (74) 



F±{t) = 1 +^/„±-r", (75) 

n=l 



00 [\/9] 

B{t) = Y.Y1 ^'^''■' • (Inr^)^ • T«, (76) 

q=0 p=a 

with r taken to be real, and where the &*^P''?^'s are the same for T < T^ and 
T > Tc, and are functions depending on the lattice. Unlike the expansions (jH]), 
and (P5)) . which converge we will see below that (|76p is expected to only be 
asymptotic. The function F±{t) is referred to, in [TT] and [33], as a "scaling 
function" . 

It is greatly instructive to compare the result (|74|) of [11] and [26] with 
the behavior of the expansions (|44| and (|45l) . derived from the form factor 
expansion. 

The term B{t) is of the same form as the logarithmic terms (f22]) already 
seen in the determinantal form of the correlations. Such terms are present term 
by term in ([5]), and are referred to, in [11], as "short distance" contributions. 
However, in the expansion (j44|) and (j45)). the x^^H'^) have singularities, at 
T -^ Tc, with powers of In \T - Tc\. Thus, term by term in dH]) and (gS]), the 



•^As far as analysis of several complex variables 1641 is concerned, the situation is even worse 
for the anisotropic Ising model, because of the codimension zero manifold S of section 3.6. 
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coefficient {l — tY''^ is not cancelled out. Consequently there must be an infinite 
number of detailed cancellations to obtain the B{t) of ([76]) from (|44|) and (|45l) . 
The most interesting question concerns the convergence of the infinite series 
in ((75|) and ([75)) . This cannot, of course, be answered by a finite series, no matter 
how long and, in fact, the numerical results for F±(t), and i^P'"^) in |2S], do not 
show divergent behavior with the number of terms which have been computed. 
However, the dense set of singularities in the x^"H^)j found by Nickel [9], [lO] 
in analyzing the form factor integrals, must have a significant influence on this 
expansion. The influence of Nickel's singularities has been analyzed in |11| with 
the conclusion that there must be asymptotic behavior in, at least, some of the 
series in ([75)) and ([7S)) . It is argued in [llj, from the results of (6.12), that 
the coefficients b^P'^^ form an asymptotic sequence for sufficiently large p. The 
behavior of b^^''^' , for p fixed and q — ?> oo, and /„± as n -^ oo, remains to 
be carried out. The further effects caused by the Landau singularities [TO] also 
remain to be studied. 

3.9 The scaling (field theory) limit at i^ = 

One of the most important properties of the Ising model at -ff = 0, as T — ^ Tc, 
is that it defines a Euclidean quantum field theory and gives a very concrete 
example of the concepts of mass and wave function renormalization. 

We concretely illustrate this for the two point correlation (cco'^a/.w)- For 
mass renormalization we define 

r = k{T; M/N) ■ R, (77) 

where R^ = M^ + iV^, and where n{T] M/N) is the inverse correlation length 
introduced in subsection (3.4). We define what we call the scaling limit, for T 
real, as 

T — > T,,, R — ^ oo, (78) 



with r fixed, and we recall from psp that k vanishes as T —>■ Tc 

k{T-M/N) ^ A^- {l~t). (79) 

For wave function renormalization we divide the correlation function by the 
factor (1 - i)!/" which vanishes at T = T^. When T < T^ this factor is the 
square of the spontaneous magnetization. The interpretation is that the spins, 
which on that lattice have the values tr = ± 1 , are regarded as having a "natural 
size" of M. in the scaling limit. Similarly, for T > Tc, the value of the disorder 
parameter M^ is interpreted as the natural size for a. 

We then define for T -^ Tc± 

G±ir) ^ lim (1 - t)-!/-* • (ao,ofTA/.jv). (80) 

scaling 

By use of the form factor expansion we see that this limit exists. In the isotropic 
case this function is rotationally invariant. In the anisotropic case it becomes 
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rotationally invariant if one uses the length variable r (fixed) 



(sm\v2E^/kBTc\^''^ ^2 (smh2E''/kBT^V'^ ^ 



M' + -^^T^T-^TT^^ • N 



-,1/2 

\amh2E^/kBTcJ ^'^ ' \smh2E'^ /keTc J " '^ ^ *> 

(81) 
This function G± (r) is expressed in terms of a Painleve equation of the third 
kind 

2 7](r/2)V-: 4 J^/2 V 

The second order equation ([5^ admits a one parameter family of solutions 
which decay exponentially, when — > oo, as 

vie) ^ 1 - — ■ Ko{2e) (84) 

TT 

where Ko{26) is the modified Bessel function. The specific value A = 1 cor- 
responds to the Ising model. The result was first announced in [56j and [65]. 
Two different proofs were originally given. The first, in |^, is based on Myers' 
work [66], on the scattering of electromagnetic radiation from a strip, and the 
second [7] is based on a direct manipulation of the exponential representation 
in the scaling limit. 

An alternative derivation of the scaled two-point function was subsequently 
obtained [S] by directly taking the scaling limit of the Painleve VI equation 
([36]). This scaling leads to representation in terms of a Painleve V function. 
The equivalence of this representation with the original Painleve III form was 
shown in |67) . 

The scaling limit as defined by (|77)) and (|75)) . which is used to define the 
scaled Green's function (|5(I1) . is defined from the massive regime where the cor- 
relation on the lattice decays exponentially. It remains to connect this regime 
with the algebraic decay of the lattice correlations, at T = Tc, given by p4)) . 
We do this by extending the strict limiting definition (IMl) with the less precise 
statement that 

{<yo^o<JM.N) - {l-ty/^- G±{kR), (85) 

and examine (|85p as T — >■ Tc, i.e. i — > 1 and k — >■ oo. Thus, we see that if, as 
r — >■ 0, we have 

G±(r) ^ AGy/\ (86) 

then (|55]) reduces to 

/A \^/^ 



16 



The exponent 1/4 in (|86l) is shown in [7] to follow, from the local behavior of 
the Painleve III equation at — >■ 0, if the exponentially decaying boundary 
condition ([5^ . at infinity, holds with A = 1. 

The constant Aq does not follow from a local property of the Painleve func- 
tion, and the limit of G±{r), as r — >• 0, will agree with (cto,oO'a/,7v), as i? — >■ oo, 
if in addition to the functional form (1551) , it can be shown that 



Ag-AI^^ = A,, (88) 

with Ac given by (j25p . This crucial identity was first proven by Tracy [68' in 
1991, by use of the scaling limit of the exponential forms of the correlation (|26p 
and (1271) . It is perhaps worth pointing out that a derivation of the constant Aq 
of (|86p has never been directly obtained from the Painleve III form 



4 The Ising model for H ^ 

The properties of the Ising model at H = 0, presented in the previous section, 
are all found by exact computations which start with the definition ([T]) of the 
nearest neighbor Ising model and are mathematically rigorous. However when 
we extend our considerations to _ff 7^ 0, this is not the case, and, with only a 
few exceptions, the results require some arguments which, while often extremely 
plausible, in fact include assumptions which remain to be proven. Nevertheless 
the work of the last 50 years has produced remarkable results which give a 
compelling scenario of the behaviour of the Ising model for H ^ 0. In this 
section will here discuss the following major contributions: 

1. Perturbation for small H for the two-point function US^ and two-body 

scattering for T > T^ [55] . 

2. Conformal field theory |34| . 

3. Integrable perturbations of conformal field theory [35], [36] . 

4. The connection of integrable perturbations with Rogers- Ramanuj an 

identities [5g]-P]. 

5. Scaling theory with irrelevant operators and nonlinear scaling 

fields [n], [26], [44], [4^, 

6. Non-integrable perturbations of conformal field theory I46j . 

4.1 Perturbation for i^ ~ 

For T < Tc the two point function has been studied perturbatively [3^ for 
small H. It was found that if in the limit i? — >■ and T — > Tc- from the low 
temperature side with the "scaled magnetic field h fixed 

TT 

h = lim —7^, (89) 

scaling T - Tc^'^/S ' ^ / 
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then the connected part of the two-point function, for h ^ and large r, is 







with 



K,j{h) = /i^/^ • Aj and aj{h) = h ■ a, (91) 



where a is a constant, and the A^'s are sohitions of 

Ji/3(A^/V3) + J-i/3(A2/V3) - 0, (92) 

where Jn{z) is the Bessel function of order n. 

This perturbation is, in fact, a "singular" perturbation which is subject 
to some interpretation. In particular one is not able to distinguish between 
r^^/^ e^*"" and Ka{K, r), where Kq{z) is the modified Bessel function. If we 
make this replacement in (IMl) . then the Fourier transform consists of a set of 
poles, and this result can be interpreted as an example of confinement of two 
"domain wall excitations" which interact by means of an infinitely weak linear 
confining potential, to produce a spectrum of "mesons" . For T > Tc the only 
effect on the leading behavior of the two-point function, for h ^ 0, is to shift 
the correlations length by a term proportional to h^ . From these computations 
a scenario is conjectured in [5?] that, as we go from T < Tc to T > Tc in a path 
in the (iJ, T) plane in the scaling limit, the Fourier transform of the two-point 
function will contain poles. These poles can be interpreted as quasi particles 
and the number of these poles will go, in a smooth fashion, from one at T > Tc-, 
H = 0, to the infinite number of poles for T < Tc and /i — > given by (|90l) . 

Much more recently [33] the two-body scattering has been studied in the 
same small h limit for T > Tc- One of the very significant results of this study 
is that, at large energies, the scattering is predominantly inelastic. 

4.2 Conformal field theory 

Conformal field theory is an entirely new approach to critical phenomana in- 
vented by Belavin, Polyakov and Zamolodchikov [33] in 1984. In this approach 
there is no lattice such as ([1]) and the theory is defined directly in the contin- 
uum. In particular, these continuum theories make no contact with the short 
distance terms, and the natural boundary, which were discussed in section 3. 

It is not our purpose here to present a survey of conformal field theory, 
which is well presented in the original paper [34] and, subsequently, in many 
places fB^, [70]. Instead, we will restrict our attention to a discussion of two 
integrable perturbations relevant to the Ising model, the M(3, 4) and the M(2, 5) 
minimal models. 

The minimal model M(3, 4) describes the critical point of the Ising model 
at T — Tc and H = 0. It has two relevant operators for the energy and magne- 
tization. The conformal dimension of energy is 1/2, which means the two-point 
function for energy is r"^, and for magnetization is 1/16, which means that the 
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two-point spin-spin correlation is r~^^'^. These results agree which the exact 
results at T = Tc of the Ising model ([1]). 

There is, however, a second conformal field theory which is relevant to the 
Ising model, namely the M(2, 5) model which is related to the Lee- Yang edge. 

In 1952 Lee and Yang [7T] proved what is one of the few results exactly 
known for the Ising model in a magnetic field, namely that, for real interaction 
energies, the zeros of the partition function of an isotropic Ising model on a 
finite lattice, all lie on the circle z = 1, where 

z = e-2ff/^-«^. (93) 

For T > Tc these zeros are all bounded away from the unit circle z = 1, and 
they pinch the real z-axis when T — > Tc- For T < Tc they fill up the unit circle 
1^1 = 1. For T > Tc the endpoint of the arc of zeros is called the Lee- Yang edge. 

The confining of partition function zeros to an arc in the complex z plane 
for real temperatures is to be contrasted with the zeros of the partition function 
in the complexplane u — q-'^^/'^bT j^j. ^^^^ H, where computer studies [72., on 
systems of sizqj up to 16 x 16 show that the zeroes for ff = are located on 
curves only for very special boundary conditions [73] , and for iJ ^ there are re- 
gions of the complex u-plane where, even for these special boundary conditions, 
the zeros do not lie on curves. 

The identification of the Lee- Yang edge as a critical point, with a continuum 
field theory interpretation, was first made by Fisher [71], and the identification 
of this field theory with the M(2, 5) minimal model was first made by Cardy [37] . 
There is only one relevant operator and the dimension is —1/5, which means 
that the two-point function of this operator is {4>o4>r) is r^^^. Because there is 
only one relevant operator, both the energy operator aj^kCj^k+i and the spin 
operator aj_k rnust have the same two-point function as the operator (j) of the 
M(2, 5) model. 

The value of the magnetic field H = i Hly, at the Lee- Yang edge, vanishes 
at T ^ Tc from above, as: 

Hly = Cly- [T-Tcf'''^. (94) 

Cardy has determined [27] that the density of zeros p{Im{H)) diverges, when 
Im{H) — >• Hly, as: 

p{Im{H)) ^ {Im{H)-HLYy^^^- (95) 

4.3 Integrable perturbations of conformal field theory 

A perturbation of a conformal field theory is called integrable if it preserves an 
infinite number of the constants of the motion of the conformal field theory which 
is being perturbed. This concept was introduced by Zamolodchikov [33], [33] at 
a conference in 1988, where he discovered that the perturbation of the M(3,4) 



*We have also performed, with I. Jensen, calculations of zeros of partition function of the 
square Ising model in a magnetic field up to size 20 X 30. 
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model, by both the energy and separately by the spin operator, have an infinite 
number of conservation laws. The M(2, 5) model was shown [38] to have an 
integrable perturbation by the operator cj) in 1990. 

4.3.1 The case of A/(3,4) 

The perturbation of the minimal model M(3,4) by the energy operator e = 
'^j,k'^j,k+i is, in fact, nothing more than saying that the original Ising model 
with nearest neighbor interactions ([T]) is integrable, which is manifestly true 
and completely not surprising. However, the discovery in [35', and [36', that 
the perturbation by the spin operator (Jj^k, which corresponds to the Ising model, 
at T = Tc in a non zero magnetic field, is integrable came as a big surprise, 
because the lattice interaction ([T]), at the critical temperature with H 7^ 0, is 
not integrable. Fortunately, this mystery was resolved in 1992 when an inte- 
grable lattice model was found [7S] which does have the critical behaviour of 
the magnetization of H^/^^, found in [33; and iHS', for the perturbed conformal 
field theory. This is the behavior of the Ising magnetization obtained, decades 
ago, by simple scaling arguments for critical exponents. 

The truly remarkable property of both, the perturbed conformal field theory 
and the lattice realization, is that they have an excitation spectrum with eight 
quasi-particles, and that the masses of these particles are proportional to the 
components of the Perron-Frobenius eigenvector of the Cartan matrix of the 
£'§ Lie algebra [35] j [33- This is completely in accordance with the scenario 
proposed from the perturbative computation of |32j . 

4.3.2 The case of M(2, 5) 

The integrable perturbation of the M(2,5) model was investigated in ^38]- In 
this case there is only a single quasi particle excitation. A lattice realization is 
regime I of the hard hexagon model j76' . 

Because there is only one relevant operator the identification of this per- 
turbation with the Lee- Yang edge, for real T ^ Tc and complex magnetic field, 
should be equivalent to the corresponding edge in the complex energy (tempera- 
ture) plan for real magnetic field H ^ 0. However, the complex energy partition 
zeros do not lie on curves the way the complex magnetic field zeros do. Thus, 
for this identification to hold, further restrictions on these complex energy zeros 
must hold in the vicinity of the edge. Consequently, the precise relation the 
energy edge has with the perturbed M(2, 5) model remains an open question. 

4.4 Roger s-Ramanuj an identities 

All conformal field theories possess a Bose-Fermi duality [41], [42]. The Bose 
form gives characters in the Roccha-Caridi form [77j by eliminating singular 
vectors from a Bosonic Fock space. For the minimal models M{p,p'), these 
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characters are 



^(P,P') — . y^ (gjipp'j+rp'-sp) _ (p'j+s)(pj+r)\ /gg\ 



where 



(qU = 11(1-'?')' (97) 

fe=i 

for 1 < r < p — 1 and I < s < p' — 1. This is a "field" representation of 
the theory, and can be thought of as being natural to characterize "short dis- 
tance" (ultraviolet) properties. An equivalent representation is given in terms 
of exclusion statistics applied to a set of Fermionic quasi-particles. This gives 
rise to a Fermionic form of the characters. These Fermionic forms give a par- 
ticle characterization of the spectrum, and can be thought of as the natural 
characterization of "long distance" (infrared) properties. The relation between 
the Bose and Fermi forms is a generalization of the famous Rogers- Ramanuj an 
identities [75]. [75]. 

The illustration of this duality for the Ising model, a.t T — Tc, H — 0, as 
the Af (3, 4) minimal model is very instructive. 

From (2.8) of [39] the Fermionic representation of character x^ 2 ' ^^r the 
spin operator of the minimal model M(3,4) which characterizes the Ising con- 
formal field theory, is 



, ^ °° „m(m-l)/2 



Xi,2 



E 



„?7i(m — 1)/2 

* • (99) 



{q)r 



(3,4) 



Similar formulas hold for the identity character xi i ; s-^id the energy character 

X2 i ■ These Fermionic forms match the particle excitations seen in the form 
factor expansion of the form factors for H = 0. There is one type of excitation, 
and the index m represents the contribution of a m-particle state. The sum 
over odd (even) m in ([55]) (respectively (|M1) ). corresponds to the odd (even) 
number of quasi-particles in the form factor expansion of the two-point function 
for T>Tc {T < Tc). The equality of dMI, and (jMl), represents the fact that, at 
T = Tc, these representations of the spectrum for T ^ Tc must give the same 
result. In the language of perturbed quantum field theory this characterization 
of the spectrum corresponds to perturbing the Af(3,4) minimal model by the 
energy operator. 

There is, however, another Fermionic representation of these same characters 
which was first conjectured in (2.18) of "401 (and proven in [43] ) 

00 00 mCp^m* 

mi=0 ms=0 *>*''"i ^^'""^ 
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where m = (rrii, • • • mg,), and where Ces is the Cartan matrix of the Lie algebra 

Es. Similar identities hold I4T] for X?f+Xi.f and Xm^^ +^1^2''^ +xi^^^- This 
representation has eight different types of excitations, which have no restriction 
on the parity of the number of allowed excitations. This is in exact correspon- 
dence with the eight particles, found by Zamolodchikov [3Bi by perturbing the 
M(3, 4) conformal field theory by the spin operator. 

For M(2, 5) the fermionic forms of the characters have only one quasi- 
particle, and read 

(2,5) _ ^g^ (2,5) _ Y^^ .^o^^ 

n— n— 

which are, in fact, the original identities of Rogers [78] and Ramanujan [79] . 

4.5 Scaling theory 

A completely different approach to the Ising model, with H j^ 0, is the scaling 
theory of the renormalization group which was developed before the discovery of 
conformal and perturbed quantum field theory. This is presented in detail the 
classic paper of Aharony and Fisher [44 , and tested extensively in JJJ and [55] . 
In [44] the conjecture is presented, for T and H real, that 

F =. -gl ■ \ngl ■ Y + g^ ■ Y±igh/\gt\''^^) + Ao{t), (102) 

where the function Ao{t) is analytic at r = 0, and where Y±{z) refer to T above 
and below Tc, and 

9t = Y. «2nW ■ ^'", ao(0) = 0, (103) 

n=0 



9h 



J2 b2n+i{r) ■ i?2"+\ (104) 



with 



n=0 



00 



a2n(T) = ^ a2„,j • r^ with ao,o = 0, 

00 

62„+i(r) = Y. ^2n+i,J ■ ^'- (105) 

j=o 

The functions gt, gu may depend on the interaction energies E{j,k) in (jlip . 
but the constant Y, and the functions Y±{z), do not depend on the interaction 
energies. The functions Y±{z) are referred to as "scaling functions". The inde- 
pendence on the constants E{j, k) is referred to as "universality" . The leading 
behavior as T ^ T^, and iJ ^ 0, of gh/\gt\^'^^^ is: 

|ao,iT|^^/« 
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(where the last hne is a definition). When H ^ the free energy is analytic 
for all real values of T, and, therefore, the functions y+{z) and Y-{z) must 
analytically connect with each other as z — > oo. These are, in fact, parts of the 
same function. This is conveniently expressed by defining a new variable 

Instead (I102p , we follow [?B] in defining the "singular" part of scaling free energy, 
for T and H real, as 

2 

F = ^- Inr^ +|if|i6/i5. $(^)^ (108) 

(where we have used the normalization conventions of |46) with the exception 
that we use — r instead of m). For H and r both real and positive, we have 

77 = l//i«/i^ (109) 

Thus, for real 77 > (T > Tc), we have 

$(77) = r,^. y+(l/r7i5/8). (110) 

For i7 > and r < (T < T^), we have 

ry = -1/h^^/^. (Ill) 

Thus, for real ry < we have: 

$(,7) = r,2-r-(l/(-r/)i5/«). (112) 

The function $(77) has, for large values of r] on the real line of [13], the 
behavior 



$(77) = ^' • £ Gfc • (-77)-i5^-/8^ ^ ^ _^^ (113) 

00 



fc=i 
and, for small values of 77: 



$(?7) = -^In?/' +^$,.•77'=. (115) 



Stt 



k=0 



For real values of 77 the function $(77) has been numerically determined |80) 
by Baxter's variational approach, based on the corner transfer matrix [ST], |82] . 
as enhanced by an improved iteration scheme [83, . 

From the form (|102l) Aharony and Fisher [44] derive an expression for the 
functions -F^ (t) of (|74l) , which is supposed to have validity, not only in the limit 
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7J — > 0, T — > 0, but also for H and r finite. Unfortunately, the conjectured 
form (|102[) fails for several reasons. First of all, the result of 04] has F+{t) = 
F-{t), whereas the analysis of [TT], and [35], show that F±{t) differ, at order 
T^, on the square, triangular and hexagonal lattice. In addition, for the square 
lattice, the term t"^ of [33] disagrees with [IT], [3SJ. Furthermore the conjecture 
(|102p does not account for the term B{t) of ([7^ . Clearly a more general 
conjecture is required. 

A more general conjecture, than that of Aharony and Fisher, utilizes the 
"irrelevant variable" of scaling theory. We follow the notation of [TT] and [55] . 
We write this more general conjecture, for the singular part of the free energy, 
as 

+.9? • Y±{g,/\g,\^'/\{g^J\g,\y^/y^}), (116) 

where, now, the coefficients a2n(i) and 62n+i, in the nonlinear scaling fields in 
(|104p , are allowed to depend on a set of variables Uj , and where 

5«. = ^c^^J^,")-^'", (117) 

are additional "nonlinear scaling fields" associated with the irrelevant operators 
of scaling theory that have dimensions yj . Further definitions of these concepts 
may be found for example in [45) . However, no explicit forms, or conjectures, 
for these multivariable formulas have ever been given. Moreover there is no pre- 
scription given to separate the effects of the "non linear scaling fields" from the 
irrelevant operators. In addition the higher powers of Inr, seen in the suscep- 
tibility, are not present in the form ()116|) . Consequently, while the form (|116p 
may be regarded as "conventional" , it is descriptive rather than computational. 
Indeed, in [TT], and [26], even though it is introduced, it is never used. 

A further difficulty with the "scaling form" (|116p is that it does not give an 
explanation for the higher powers of In \t\ which occur in B{t) of (|22|) . At this 
point the literature is slightly ambiguous. In [TT], in footnote 12, it is stated that 
the fact that there is a "resonance between the identity, and the energy" , will 
result in higher powers of In r (much as integer differences of exponents lead to 
powers of logarithms in Fuchsian differential equations). However, in footnote 5 
of [in], it is said that these powers of In |t| have "not yet been interpreted within 
the context of scaling theory" . There are clearly things left to be understood. 

To see what is needed, we generalize our point of view, from the nearest 
neighbor Ising model ([T]), to the much more general case ([TT]). which allows for 
many further neighbor interactions. This more general model is not integrable. 
Now the internal energy p^ will include the spin correlations of all spins which 
are connected by non zero interaction energies. Thus, at least for small values 
of the non planar bonds, it is entirely reasonable to expect that, at T —>■ Tc, the 
singularity will not be a pure logarithm, as is the case for the nearest neighbor 
case, but will involve many higher powers of In^ |t|, each of which are multiplied 
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by some suitably high power of r. We might speculate, from ([2^ . that this 
power will be r^ , but there is no further argument for this, except that it is 
the result found in (22). Moreover, even in the case of the nearest neighbor 
interactions ([T|), these terms, with higher powers of In |r|, should be present. By 
including such terms we will be able to reproduce the "short range" term B{t) 
in the susceptibility expansion ((76)) . 

But the truly serious problem in obtaining a general form of the free energy 
of the Ising model in a magnetic field, as an expansion for "small H and small r" , 
is that the natural boundary, discussed in section 3 for the susceptibility, forces 
the expansion ([7i)) to be asymptotic. Therefore, even though the susceptibility 
is well defined for all values of r, the expansion ([74| cannot be taken as the 
definition of the susceptibility. This is exactly the same problem which afflicts 
perturbation expansions in the Feynman diagram expansion of quantum field 
theories, such as Quantum Electrodynamics [3^, [5T| . 

4.6 Non-integrable perturbations for H ^ 0, T ^ 

The most ambitious program of studying the Ising model, for H ^ 0, is to allow 
both T 7^ Tc, and H ^ 0, at the same time. This will allow the passage, from 
T < Tc, if = to T > Tc, if = 0, on a path in the real (ii,r) plane, as is 
contemplated in the scenario of [31]. However this more general two variables 
perturbation does not satisfy the criteria of [3S], [33] needed for an integrable 
perturbation. Furthermore, in this more general case, we cannot treat the region 
of if ~ 0, and T ^ Tc in isolation from the Lee- Yang edge for complex H. 
Thus, we need to understand how the Af(3,4) conformal field theory fiows to 
the Af (2, 5) theory under the infiuence of the nonintegrable perturbation. 

To make this precise it is necessary to define a scaling limit in the complex 
plane. This is very delicate because it says that we are describing the singu- 
lar behaviour of a function of two complex variables in terms of one complex 
variable. Such a reduction will require very special circumstances to be valid. 

Some of the consequences of the existence of this limit have for the free energy 
are extensively discussed in the 2003 paper of Fonseca and Zamolodchikov |46] , 
where a scenario is presented which incorporates several assumptions about 
analytical continuation in T and if, as two complex variables. We note, however, 
that the question of the relation of the natural boundary in the complex T plane 
to the analyticity properties in the complex H plane near the Lee- Yang edge 
remains to be investigated. 

5 Why the Ising model is important 

It is very natural to extend perturbed conformal field theory of the Ising model 
from the integrable cases of [3S], [33] j where either H = 0, T ^ T^, and 
T = Tc, if ^ 0, to the general case T ^ T^, H ^ 0. This has been done 
in [?B]. However, in spite of the impressive results of [5D] for the free energy. 
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which are in correspondence with the picture of [IS], there are several very 
interesting questions which remain to be addressed. 

Perhaps the foremost of these questions is what may be a significant differ- 
ence between integrable and non-integrable systems. We saw, in section 4.4, 
that the characters of the conformal field theory of the M(3, 4) minimal model 
have three different representations which are in one to one agreement with the 
spectrum of excitations, in the massive case, away from the critical point. We 
can thus think of the integrable perturbations as being precisely tuned to these 
three different bases, in the same way as degenerate perturbation theory picks 
out a distinguished basis. 

However, in the non-integrable cases, the spectrum of excitations contains a 
number of particles which depends on both T and H, and is not just either one 
or eight. The implication of this mismatch between the variable number of exci- 
tations, seen in the massive model (in what can be considered the infrared part 
of the spectrum), and the existence of only three particle- like representations 
of the character in the conformal field theory (which can be considered as the 
ultraviolet part of the spectrum) , is that the smooth match, found in the Ising 
model at H — (which was shown in section 3.8 with the precise demonstration 
by Tracy ^68; of the equality ([55]) ). may not hold in the general non-integrable 
case. 

The discussion of the potential disagreement of the short distance behavior 
of the scaled two-point correlation function for general values of T and H of 
the Ising model, with the result expected from conformal field theory, raises the 
suggestion that our understanding of scaling theory may not be complete. 

Scaling theory is concerned with the relation of two length scales: the scale 
of the lattice length on which the theory is defined, and the correlation length 
which is observed in the correlation functions. In scaling theory, as used in 
statistical mechanics, and in renormalization of quantum field theory, there is 
a smooth match between these two scales. But this is drastically different from 
what occurs in systems, such as fluid mechanics, where a common piece of folk 
wisdom of perturbation theory [M] is that "divergence should be expected when 
the solution depends on two independent length scales" . This would also seem 
to be related to the fact, found in llj and [55], discussed in section 3.8, that 
the susceptibility does not have a convergent expansion about T — Tc- It would 
also seem to be in agreement with the fact that the natural boundary of the 
susceptibility, presented in section 3.7, does not naturally fit into conventional 
scaling theory. Thus it may be the case that the smooth matching of long 
and short distance expansions of the Ising model, at if = 0, is caused by the 
integrability of the model, and may not hold in the general non-integrable case. 

For these reasons it can certainly be said that, even though the Ising model 
is the best understood system in statistical mechanics, there are still many 
puzzling questions to be investigated, questions which have an importance well 
beyond the narrow range of just this one model. In fact, it can be argued that 
many of the nagging questions concerning our understanding of quantum field 
theory [30], [31] are related to the puzzles of the Ising model. Hopefully these 
questions will inspire future research. 
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